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We consider the perturbative description of saturation 
based on the nonlinear QCD evolution equation of Balit¬ 
sky and Kovchegov (BK). Although the nonlinear corrections 
lead to saturation of the scattering amplitude locally in im¬ 
pact parameter space, we show that they do not unitarize 
the total cross section. The total cross section for the scat¬ 
tering of a strongly interacting probe on a hadronic target 
is found to grow exponentially with rapidity t = ln{s/so), 
a oc exp{ et} where e is a number of order unity. The 
origin of this violation of unitarity is the presence of long 
range Coulomb fields away from the saturation region. The 
growth of these fields with rapidity is not tempered by the 
nonlinearity of the BK equation. 

Understanding the growth of total scattering cross sec¬ 
tions with energy ^/s is a longstanding problem. The uni¬ 
tarity, or Froissart bound states that the total inelastic 
cross section for the scattering of a hadronic projectile 
on a hadronic target can not grow faster than 

CT < 7rd^ ln^(s/so) , (1) 


where U (x) is the eikonal scattering amplitude of a fun¬ 
damental probe on the target characterized by some dis¬ 
tribution of gluon fields A^. In the gauge used in 


U{x) = P exp 


dx"T“A+(x) 


( 3 ) 


The averaging in (i is taken over the ensemble charac¬ 
terizing the target. In the large Nc limit, (^) simplifies to 
a closed equation for the scattering probability N{x,y) 
of a colour singlet dipole with charges at points x, y 

N{x,y) = ^Tv{l-U\x)U{v)). (4) 


This equation was independently derived by 
Kovchegov [Q in the colour dipole approach of |Q| 


4iV(x,!,) 

[fV(x, 


z) + N{y, z) - N{x, y) - N{x, z)N{z, y)] . 


( 5 ) 


where d is some typical hadronic scale and t = ln(s/so) 
the rapidity. While QCD, the theory of hadronic interac¬ 
tions, is a unitary theory and therefore satisfies this uni¬ 
tarity bound, there is no guarantee that perturbative cal¬ 
culations preserve this property. In fact, the linear per¬ 
turbative evolution equation due to Balitsky, Fadin, Ku- 
raeev and Lipatov (BFKL) implies an exponential growth 
of a with t, thus violating unitarity. 

Following the pioneering works of ||^,|| , there has been 
recent progress in high energy hadronic scattering in a 
derivation of a nonlinear evolution equation ti which 
tames the BFKL-type growth. These equations resum 
the nonlinear corrections to the QCD evolution with ra¬ 
pidity (energy) to all orders in partonic density and to 
first order in the QCD coupling. While previous stud¬ 
ies of these equations assumed translational invariance 
in the impact parameter plane, here we explore for the 
first time their impact parameter dependence. For the 
equations first derived by Balitsky 1^, we show that the 
total cross section does not unitarize but grows exponen¬ 
tially with t. 

The first in the hierarchy of nonlinear BK evolution 
equations which govern the evolution of correlation func¬ 
tions of the gluon helds in the target with t, reads 


At, a - uH.)u(y)) = 


{x-y^ 


(x - z)"^{y - z)2 
{N,Ti[U\x)U{y)] - Ty[UHx)U{z)]Tr[U\z)U{y)]) , (2) 


Weigert suceeded to reformulate Balitsky’s hierarchy 
in terms of a nonlinear stochastic process. 


dU (a;) 
dt 


9 U{x)iT 

io-s f 

27T^J 


f d^z {x — 

(x — 
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(x — z)2 


z)i 

Z)2 


l-U'^(x)U(z) 



TtIT‘‘U^x)U(z)]. (6) 


Here U (x) and 17(x) are the unitary matrices (^) in the 
fundamental and adjoint representations, respectively. 
The noise is characterised by Gaussian local correlations 

(e“(t', z")) = - t")5{z' - z") . (7) 


This simple Langevin equation gives rise to an infinite 
number of equations for correlators of U which coincide 
with those derived in . 

Nonlinear evolution equations were also derived in the 
Wilson renormalization group approach B. While re¬ 
sults from this approach coincide with (^,^ in a certain 
limit, the question whether they are generally equivalent 
or incorporate different physics is still open. Here, we 
focus entirely on the BK equations (||,^,||). 

From the first numerical and analytical 00 

studies of the BK eqs. (^,^ the following consistent pic¬ 
ture emerges: Suppose one starts the evolution from the 
initial condition of small target helds (or N{x, y) <C 1 for 
all x,y). Then initially the evolution follows the BFKL 
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equation, since the nonlinear term in is negligible. As 
the scattering probability approaches unity, the nonlin¬ 
ear term kicks in and eventually the growth stops as the 
RHS of (||) vanishes for N(x,y) = 1. The larger dipoles 
[(x — » 1/Qs(0] saturate earlier, the smaller dipoles 

follow at later ’’time” t. These features are contained in 
the simple parametrization |T^ 


N{x,y) = 1 — exp 


{x-yfQl{t) 


( 8 ) 






(11) 


This equation neglects contributions to the derivative of 
U that come from gluons originating from the sources 
outside the black region. Those contributions would en¬ 
hance the growth of U, and so by omitting them we un¬ 
derestimate the rate of growth of the radius of the black 
region. The formal solution of eq. & is 


The saturation momentum Qg (t) is a growing function of 
rapidity. Its exact dependence on rapidity is not known, 
but both, the numerical results Q and simple theoreti¬ 
cal estimates |]14[|12f| are consistent with the exponential 
growth of the form 


Qsit) = A exp 


asct 


( 9 ) 


1-U{x,t) = 




/ 

/to J\z\<R{t) 


2 {x z), 


{x — z)- 


U^)- (12) 


Squaring it and averaging over the noise term gives 


{{i-u{x,t)r) = 


sNc 


dr 




'to 


<\z\<R{r) 


{x — z)- 


• ( 13 ) 


with c of order unity. This physical picture has been 
anticipated several years ago in Q. 

While the BFKL equation leads to an unphysical expo¬ 
nential growth of the scattering probability N{x, y) with 
t, the nonlinearities of eqs. (||,||) tame this growth such 
that N{x,y) < 1, as required for a probability. This is 
commonly refered to as “unitarization”. However, the 
“saturation” of the scattering probability at fixed im¬ 
pact parameter does not insure that the total scattering 
cross section is unitary. We hence refer to the above phe¬ 
nomenon more accurately as ” saturation”. 

To calculate the total inelastic cross section one has 
to integrate the scattering probability over the impact 
parameter. Thus in the saturation regime 

(T = 7ri?^(t), (10) 

where R{t) is the size of the region in the transverse 
plain for which the scattering probability for hadronic 
size ’’dipoles” is unity. To satisfy the Froissart bound 
the radius R{t) should grow at most linearly with t. We 
now present two simple calculations which establish that 
within the BK evolution the growth of the radius with 
rapidity is exponential. 

First consider the Langevin equation (i- Assume that 
initially, at rapidity /q the target is black within some ra¬ 
dius Rq. This means that for \z\ < Rq the matrix U{z) 
fluctuates very strongly so that it covers the whole group 
space. We concentrate on a point x which is initially out¬ 
side of this black region. The matrix U(x) then is close 
to unity. Thus there is no correlation between U(x) and 
U(z), and the second term on the right hand side of (|^) 
can be set to zero. This is the random phase approxima¬ 
tion introduced in and used later in . As the target 
field ensemble evolves in rapidity, the radius of the black 
region grows. As long as the point x stays outside the 
black region we can approximate the Langevin equation 
by (we drop colour indices which are inessential to our 
argument) 


As long as x is outside the black region we can approxi¬ 
mate the integral on the right hand side by 




RHr) 


J\z\<Rir) {X-Zf X^ 

and eq.(p^ becomes 

{{1-U{x,t)f) = f dTR^ir). 


(14) 


(15) 


Now as the black region grows, eventually it reaches the 
point X. At this rapidity the matrix U(x) will start fluc¬ 
tuating with the amplitude of order one. Thus when 
R(t) = |a;|, the left hand side of eq. ( |T5| ) becomes a num¬ 
ber of order one, which we call 1/e. We thus have an 
approximate equation for R(t) 


-R^(t) = f dTR^^r). (16) 

e Jto 

At large rapidities therefore the radius of the black region 
is exponentially large 


m 


R(to) exp 


asNc 

2tt 


e{t - to) 


(17) 


This is our main result. 

We note that while the approximations leading to eq. 
(0) cease to be valid when the point x is on the boundary 
of the black region, this does not affect our main conclu¬ 
sion. First, eq. & is an underestimate of the integral, 
thus underestimating the growth of R. Second, when x 
is on the boundary of the black region and z in the black 
region, although the factors {l — U{x)W{z)) and U{x) in 
eq. (||) are not strictly unity, they are still of order one for 
allmost all points z. Thus, although we can not deter¬ 
mine the exact numerical value of e, the functional form 
of the solution as well as its parametric dependence is 
given correctly by eq. (EH)- 
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Note that eq. refers to the evolution of matrices 
U (x). This can be thought of as evolution of the scatter¬ 
ing amplitude of a coloured probe. The preceding deriva¬ 
tion thus refers to the growth with rapidity of the cross 
section for scattering of a coloured probe. In a confin¬ 
ing theory this is not a physical quantity. However the 
physics of BK equation does not incorporate effects of 
confinement, and therefore the cross section for a colour¬ 
less dipole within the BK framework must grow in the 
same way. To establish this point, and to make more 
explicit the relation between e and the BFKL dynamics, 
we now present an alternative derivation of eq. (p^. 

To this end we consider the BK evolution as the evo¬ 
lution of the projectile @]. Suppose at the initial energy 
the projectile is a colour dipole of size Xq. It scatters on 
a hadronic target of some size i?target- As is explicit in 
[Q|, as the energy is increased the projectile wave func¬ 
tion evolves according to the BFKL equation. Thus at 
rapidity t the density of dipoles of size x at transverse 
distance r from the original dipole is given by the BFKL 
expression (see for example |p^): 

32 In 

( 18 ) 


Lut — In 


16r^ 

XqX 


MU 

XqX 

aH 


1 ^ 

and as such is valid only for In < agt. This condition 
is not satisfied by eq. However, even beyond the 

saddle point approximation the density has the form 


n{xo,x,r,t) = 


, 16r^ 

■ in-exp 

XqX 


In 

a.t 


( 21 ) 


The relevant condition is K = 0. Thus, while our calcula¬ 
tion does not specify the numerical value of e, the correct 
solution parametrically is the same as di)- 

In the target rest frame, this violation of unitarity by 
the BK evolution can be understood as follows: Start 
with a single dipole scattering on the hadronic target 
of transverse size i?target- With increasing energy the 
projectile dipole emits additional dipoles strictly accord¬ 
ing to the BFKL evolution. The density as well as the 
transverse size of the projectile state thus grows. The 
increase in density leads to increasing importance of mul¬ 
tiple scatterings which are properly accounted for in the 
BK derivation. This ensures that the scattering proba¬ 
bility saturates locally. In the saturation regime, as long 
as the size of the projectile state R{t) is smaller than the 
target size i?target, the cross section grows essentially due 
to surface effects. 


a = ttR 


2 

target 


“h 27r^target^0 


as Nc 
27T 


-et 


( 22 ) 


with oj = 41n2Nj,Q!s/7r and = lA({3)Ncas/TT. When 
the density of dipoles at a given impact parameter is 
greater than one, multiple scatterings become important. 
Thus the scattering probability is not proportional to n, 
but is an infinite series containing all multiple scattering 
terms |4j. 

For our argument, it is only important that once the 
density of dipoles at some impact parameter r becomes 
larger than some hxed critical number, the scattering am¬ 
plitude at this impact parameter saturates. The exact 
value of this number depends on the target, but impor¬ 
tantly it does not depend on rapidity. Thus the total 
cross section is given by the square of the largest impact 
parameter at which the dipole density in the projectile 
wave function is of order unity. In order to estimate this 


directly from eq. (18), we must choose the dipole size x 
in (^^ to be the smallest size which is saturated on the 
target at initial rapidity. Within the ansatz of eq. (||) 
this would he X = We then find 


R\t) = 


1 


Xo 


16 Qs{h) 


.asNc ,, 
exp[- et\ 


e = 7C(3) -1-k v'l-k81n2/7C(3) 


(19) 

( 20 ) 


Thus, also for a coulour singlet projectile, we arrive again 
at the exponential growth of the cross section. 

The exact value of e given in eq. (|^) should not be 
taken too seriously. The explicit form of the dipole den¬ 
sity eq. (^ was derived by a saddle point integration, 


R 

Thus as long as agtt < In—the cross section is 
practically geometrical. However once the energy is high 
enough so that the projectile size is larger than that of 
the target, the total cross section is determined by the 
former and grows exponentially with the logarithm of 
energy according to eq. (p^. 

This also illustrates that the applicability of the BK 
evolution crucially depends on the nature of the target. 
If the target is thick enough, so that the multiple scatter¬ 
ings become important before the growth of the projectile 
radius does, and if the target is wide enough, so that sat¬ 
uration occurs before the projectile radius swells beyond 
that of the target, then there is an intermediate regime in 
which the inelastic cross section remains practically con¬ 
stant and equal to . Then BK applies. However, 

if the target is a nucleon, neither one of these conditions 
is satisfied. Thus the tainted infrared behaviour of the 
BFKL evolution of the projectile will show up right away 
and will invalidate the application of the BK equation. 

In order to discuss the violation of unitary from the 
point of view of the evolution of target helds, we now 
go back to the stochastic process (|)- The RHS of (|) 
describes the total Coulomb (Weizsacker-Williams) field 
at point X due to the colour charge sources at points 
z. Since the noise is stochastic, the colour sources are 
completely uncorrelated both in the transverse plain and 
in rapidity. For this random source, the square of the 
total colour charge is proportional to the area, and this 
is precisely the factor R^ in eq. dil)- The incoming 
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dipole thus scatters on the Coulomb field created by the 
large incoherent colour charge. Because the Coulomb 
field is long range, the whole bulk of the region populated 
by the sources contributes to the evolution and leads to 
rapid growth of R. If the field created by the sources was 
screened by some mass, the evolution would be unitary. 
To illustrate this point, we substitute the Coulomb field 
(x — z)ij(x — z)^ in (eq.^ by an exponentially decaying 
field TO exp {—mix — z\}. It is straightforward to perform 
now the same analysis as before. Eq. (14) is replaced by 


(fzm^ exp{— to|x — z\} = exp{— to|x — R\} . (23) 

\z\<R{r} 

This leads to the substitution R^ —> explmi?} in all sub¬ 
sequent equations with the end result that 

e 


R{t) = as—t, 

TO 


(24) 


which in fact saturates the Froissart bound. Thus the 
reason for the violation of unitarity is that the evolution 
is driven by the emission of the long range Coulomb field 
from a large number of incoherent colour sources in the 
target. 

Cutting off the Coulomb field is not the only possibility 
to cure this problem. Another option is that the sources 
of the colour charge in the high density regime cease to 
be incoherent. If they have correlations ensuring that 
the total colour charge in a region of fixed size L is zero, 
then the incoming dipole would feel the Coulomb field 
only within the fixed distance L from the black region. 
Thus the new charges produced by the evolution would 
only ” split off’ the edges of the black region rather than 
from its bulk. This scenario is equivalent to exponential 
decay of the field, and will lead to a unitary evolution. In 
a confining theory like QCD, it is likely to be material¬ 
ized. We note that the desirability of such colour charge 
correlations was stressed in a somewhat different context 


Although such charge correlations do not arise in the 
BK evolution, it is not a priori clear that they are not 
present in a more complete semiperturbative framework 
which still does not take into account the physics of con¬ 
finement at low energies. In fact, the BK framework is 
incomplete inasmuch as it takes the evolution of the pro¬ 
jectile wave function to be pure BFKL. One expects that 
once the density of gluons in this wave function becomes 
large, interactions should lead to saturation effects on the 
wave function level, i.e., the density of the dipoles should 
grow slower than eq. (|T^). Such corrections should still 
be semiperturbative, in the sense that they are present 
at small as- However, for scattering on ’’small” targets, 
they will become important at the same energy as the 
multiple scattering terms resummed in eq. (^,||). These 
wave function saturation effects may lead to charge cor¬ 
relations of the type necessary to unitarize the total cross 
section. 


We finally note that the unitarity bound for Deeply 
Inelastic Scattering is different. In this case, the projec¬ 
tile is a virtual photon without fixed hadronic size. For 
transverse polarization, its perturbatively known wave 
function is $^(r) oc aem:^ for ^ Q~^■ In such a pro¬ 
jectile not all dipoles saturate at the same energy. The 
main contribution to the scattering probability comes 
from dipoles of size r > which are saturated. At 

high energy (i.e. for < Q“^), 

f d^r4>^(r) oc Q!em InQs/Q . (25) 

With the exponential dependence (||) of Qs on rapidity 
this translates into N{'y*) ex aemOts Ins/so, and therefore 

ctdis oc aemOisT^R^[t)t. ( 26 ) 

The DIS cross section has an extra power of t relative to 
the cross section of a purely hadronic process. This extra 
power of t is consistent with the numerical results of . 
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